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We study the internal structure of a non-Abelian vortex in color superconductivity with respect 
to quark degrees of freedom. Stable non-Abelian vortices appear in the Color-Flavor-Locked phase 
whose symmetry S?7(3)c+l+r is further broken to S'c/(2) c+ l+r ® ?7(1) c +l+r at the vortex cores. 
Microscopic structure of vortices at scales shorter than the coherence length can be analyzed by the 
Bogoliubov-de Gennes (B-dG) equation (rather than the Ginzburg-Landau equation). We obtain 
quark spectra from the B-dG equation by treating the diquark gap having the vortex configuration 
as a background field. We find that there are massless modes (zero modes) well-localized around 
a vortex, in the triplet and singlet states of the unbroken symmetry S(7(2) c +l+r ® c/(1) c +l+r- 
The velocities Vi of the massless modes (i — t,s for triplet and singlet) change at finite chemical 
potential /j, ^ 0, and decrease as /i becomes large. Therefore, low energy excitations in the vicinity 
of the vortices are effectively described by 1+1 dimensional massless fermions whose velocities are 
reduced Vi < 1. 



I. INTRODUCTION 



The structure of vortices is crucial in understanding the dynamics of the phases in a broad subject of physics from 
condensed-matter physics to high-energy physics and astrophysics. Macroscopic dynamics of vortices is described by 
a spatially dependent order parameter for a certain broken symmetry. It is well described by the Ginzburg-Landau 
theory, which has conducted great progress in studies of vortices in superconductors at the length scales longer than 
the coherence length and the penetration depth. On the other hand, in order to understand short range structure 
of vortices such as the vortex core, one needs a microscopic theory in which vortices are described by the fermionic 
degrees of freedom. Such a self-consistent description including both the order parameter and fermions is known as 
the Bogoliubov-de Gennes (B-dG) equation The first investigation of the vortex structure by solving the B-dG 
equation was done in Ref. [2| . Fermions are trapped inside the vortex core if their energies are less than the gap of the 
ground state while they are scattering states if their energies are larger than it. Since this first study treated the order 
parameter (the superconducting gap) as a background field (and thus it is not self-consistent), successive researches 
were devoted towards finding the self-consistent solutions to the B-dG equation d, 0] . The complete self-consistent 
description was achieved by considering both quasi-bound and scattering fermions in the vortex y| Q . These studies 
predicted an enhanced local density of fermion states at the Fermi level around the vortex core, which is experimentally 
observed in various materials 0, Q ■ Recently the analysis of vortices in the B-dG equation has been also applied to 
the BEC-BCS crossover in fermionic cold atom systems [H, [l(| . 

Among all fermion modes associated with vortices, zero modes (gapless fermions) play an especially important 
role in various subjects of physics. In non-relativistic theories including all condensed matter systems, fermion "zero 
modes" trapped in a vortex core naturally appear in the semi-classical approximation, but they are actually not 
necessarily exactly zero-energy states because they in general acquire small gaps ~ A 2 /Ep (called "minigaps") that 
are small enough compared with the gap A (Tlj . For instance, a conventional s-wave superconductor is such a case 
Q. In chiral p-wave superconductors Major ana fermions trapped inside a vortex core have exactly zero energy, which 
attain non-Abelian statistics [HI, [l3[ and give a candidate of quantum computations. Recently it has been found that 
even a conventional s-wave superconductor allows fermion zero modes in a vortex core when it is surrounded by a 
topological insulator [ijj]. Unlike non-relativistic theories, fermion zero modes in relativistic theories have ex actly 
zero energy, and the number of zero modes is determined by the index theorem to be 2n for n winding vortices [16j . 
Fermion zero modes are also important in the context of cosmic strings [l^ |. 

The study of vortices has been extended to color superconductivity in dense quark matter, where the attractive force 
between quarks induces a rich structure of symmetry breaking patterns [l8l |. At extremely high densities, a novel 
phase called the color-flavor-locked (CFL) phase will be realized, where up, down and strange quarks participate 
in the Cooper pairing with symmetry breaking pattern SU(3) C ® SU(3)l <8> SC(3)r — > SU (3) c +l+r [la420l ]. The 
color-flavor structure of the gap in the CFL phase is given by A ai cx e a ^e iik {ijj^ijj^ k ), where i, j, k are flavor 
indices, and a, /J, 7 are color indices. It was argued that superfluid vortices would exist as a result of breaking of 
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the global U(1)b symmetry of the baryon number [2l| - [23l ] . Color flux tubes have been suggested [H, [HI but they 
are topologically and dynamically unstable. The true ground states of vortices in the CFL phase are non-Abelian 
vortices found by Balachandran, Digal and Matsuura [2J] which are topologically stable superfluid vortices carrying 
color magnetic fluxes inside their core. It is a 1/3 quantized vortex, namely it has 1/3 U(1)b winding (circulation) 
and 1/9 tension of those of a global U(1)b vortex }2l| - |23j . and it also has 1/3 amount of the color flux of that of 
a topologically and dynamically unstable color flux tube [H], IH|. The detailed gap profile function and the size 
of color flux has been calculated in Ref. [Ill in the Ginzburg-Landau approach. One characteristic property of non- 
Abelian vortices is that the color-flavor-locked symmetry 5 , ?7(3) c+ l+r of the ground state CFL phase is spontaneously 
broken down to its subgroup 5J7(2) c +l+r <8> U(1) c +l+r in the presence of a non-Abelian vortex. According to 
this symmetry breaking, Nambu-Goldstone zero modes appear around the vortex, which parametrize a coset space 
CP 2 ~ 5J7(3) c+ l+r/[5 , /7(2) c+l+R (8) J7(1) c+ l+r] [H,|23|- Therefore, the vortex solution allows a continuous family of 
solutions with degenerate energy, corresponding to CP 2 . These modes are called the internal orientational zero modes 
and each point of CP 2 corresponds to the color flux which the non-Abelian vortex carries. Since this breaking occurs 
in the vicinities of the vortex, these modes are localized around it, and one can construct a 1+1 dimensional effective 
theory on the vortex world-sheet by integrating the original action over the vortex codimensions [28j . The interaction 
between two non-Abelian vortices at large distance is mediated by fluctuations of J7(1)b Nambu-Goldstone modes, 
yielding the universal repulsion between them irrespective to their CP 2 orientations in the internal space [U Hl| ■ 
This predicts that each global f (1)b vortex is dynamically unstable to decay into three non-Abelian vortices which 
have different color fluxes with the total flux cancelled out. Furthermore, it also predicts a non-Abelian vortex lattice 
at least when the lattice spacing is much larger than the vortex core size [13] ■ Such a vortex lattice is expected to be 
realized in a CFL quark matter core of a rapidly rotating neutron (or compact) star pjol l3l|. 

However, these studies are all based on a macroscopic theory, namely the Ginzburg-Landau model which is valid 
only at the scale larger than the penetration depth and coherence length. In order to understand the whole structure 
of the non-Abelian vortices including the region inside their cores, we have to consider fermion dynamics from the 
B-dG equation. The most important fermion modes are the zero modes as mentioned which exist exactly at the Fermi 
energy. The purpose of the present paper is to study the fermion zero modes in a non-Abelian vortex, based on the 
B-dG equation. 

This paper is organized as follows. In Sec. |H]we construct the fermion zero modes around a vortex background in 
the cases of a single flavor and the CFL phase. We find four fermion zero modes for a non-Abelian vortex in the CFL 
phase, which belong to the singlet and triplet representations of the unbroken symmetry S i J7(2) c +l+r- We find an 
interesting layer structure of the zero modes in the core of the vortices; the core size of the singlet zero mode is the 
half of the one of the triplet zero modes. In Sec. IIIII we construct the 1+1 dimensional effective field theory of the 
fermion zero modes. We find that the velocity of the zero modes propagating along the vortex line depend on the 
chemical potential, and that the velocity of singlet is about quarter of that of the triplet at high baryon density. The 
last section is devoted to summary of the results and discussion of the future problems. 

II. BOGOLIUBOV-DE GENNES EQUATION IN COLOR SUPERCONDUCTIVITY 

A. Single flavor 

As the simplest example which shares the essential physics with the CFL case, we first discuss a relativistic fermion 
system with a single flavor that exhibits superconductivity and allows for Abelian vortices. We treat a single (or, 
an isolated) vortex whose core is located on an infinite straight line along the z axis, and assume that the winding 
number (vorticity) is one [Hf. As long as we treat the gap profile function A(r) as a background field, the description 
in the present subsection is very similar to that of the vortex- fermion system discussed in Refs. [15l4l7| . However, we 
consider here a fermionic matter at finite densities (otherwise color superconductivity does not take place) while the 
previous studies 

tuna were formulated only in the vacuum. As we will see later, the effects of finite densities are 
quite important for low energy excitations of the fermion modes in a vortex. 

The Bogoliubov-de Gennes (B-dG) equation is useful when one considers inhomogeneous superconductivity in 
the presence of normal-super boundaries and impurities. Among such examples is the vortex in which the gap is 
locally reduced compared to the value in the homogeneous ground state. While the Ginzburg-Landau equation is 
usually used to find a vortex profile, it is valid only for spatial variations at scales longer than the coherent length. 
On the other hand, since the B-dG equation is represented by microscopic (fermion) degrees of freedom, it can in 
principle describe the internal structure of a vortex at shorter length scales. For the present case, it is expressed in 
the Nambu-Gor'kov (particle-hole) representation as 



(i) 
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with the Hamiltonian density in the mean-field approximation 

%=( -*7o7 ' V - fi A(x)7o7 5 \ ^ 
\ ~ A *(a ; )7o75 -«7o7'V + Ai/ 

Here, we consider a relativistic massless fermion and the energy is evaluated with respect to the Fermi surface (fj, is the 
chemical potential of a fermion). The gap is supposed to take a vortex profile : A(x) = |A(r)| e with |A(r = 0)| = 
and |A(r = oo)| = |A| (where r — \/ x 2 + y 2 and 8 is the angle around the vortex). 

As the first attempt to this problem, we treat the gap function as a background field as was done in Ref. @). 
However, strictly speaking, the B-dG equation ([TJ should be solved self-consistently because the gap profile function 
is expressed in terms of fermion degrees of freedom. Such a self-consistent description has been discussed in a simpler 
configuration like a kink, and recently it was achieved even for vortices in a standard type-II superconductor [5[ . We 
leave such a description in the color superconductivity for future problems. 

With the cylindrical shape of the vortex profile, one can write the eigenstates of Eq. ([1]) as ^** m (r, 8, z) with 

^ m (r,0,z) = ^ ± , m (r,e)e ik ^, (3) 

which has the momentum k z in the z direction, and chirality ± corresponding to right and left. The Nambu-Gor'kov 
structure of 1 $±, m (r,d) is expressed as 

-(£:[$>)■ ,4 > 

where particle <p± y7n and hole r] Tv m-i components are eigenstates of the z component of the total spin operator J z 
with their eigenvalues m + 1/2 and (m — 1) + 1/2, respectively [l5j]. We note that the chiralities of a particle and a 
hole (defined by the eigenvalues of 75) are ± and =F for right and left modes, respectively. 

Once the concrete form of the vortex profile |A(r)| is given, one can explicitly solve the B-dG equation ([T]). However, 
we are rather interested in low energy excitations which may be independent of the precise form of the profile. In fact, 
even without knowing the precise form of the profile, one can address the questions whether there exist zero modes 
or not, and if they exist, what are the properties of the zero modes. Note that these are highly nontrivial questions: 
in a homogeneous superconductor, even though the original fermion is massless, there is no excitation around the 
Fermi surface below the gap energy. Of course, the B-dG equation reduces to the standard Bogoliubov equation for 
a spatially constant gap, and there is no massless excitation in that case. On the other hand, since it is known that 
the vortex- fermion system in the vacuum allows for zero modes which are trapped in a vortex [15l |17| , and that its 
existence is guaranteed from the topological considerations [Iff, one naturally expects that there exist zero-mode 
solutions even at finite densities. This is indeed the case. The B-dG equation at finite densities ([T]) for a generic 
profile function |A(r)| allows for zero-mode solutions, namely, solutions with the £ = eigenvalue. The zero modes 
are realized for m = and k z = 0, and the explicit forms can be written as, for the right mode (+ for a particle, — 
for a hole) 

-Cr**™" . (5) 

n-Mn* (6) 
and for the left mode (— for a particle, + for a hole) 

= Ce-JT '«-'>l d - ( ) , (8) 

where C and C are normalization constants, and J n (x) is the Bessel function. We have represented the solutions in 
the Weyl (2-component) spinors. In obtaining the solutions above, we have imposed finiteness at r = and in the 
limit r — > 00. When // = 0, these solutions recover the zero- mode solutions in the vacuum found in Ref. [l5j], as they 
should. Notice that the zero- mode solutions above satisfy the following "Majorana-like" condition [36[ (n = ±1): 

* = "=( 7 ° 2 7 o) , (9) 
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which leads to U^WA = —H* and thus physically implies the equivalence between a particle and a hole. One can 
explicitly check that if the solution satisfies this condition, it indeed leads to m = 0, k z = and £ = 0. 

In the limit of small size of the vortex, in other words, at a distance much away from the vortex, |A(r)| can be 
regarded as a constant |A| which is equal to the gap in the bulk space. Then, for example, the right mode is given by 
a simpler form 

This asymptotic behavior clearly shows that the zero-mode solutions are well-localized around the vortex line. 



B. Color-Flavor-Locked phase 

1. B-dG equation with a non-Abelian vortex 



Now we discuss the fermion structure of the non-Abelian vortex which appears in the CFL phase. The gap 
configuration of the non-Abelian vortex in the color-flavor representation A al oc e a/37 e' jfc (V>^ T C'75'0 7fc ) (C 
is the charge conjugation, a, /J, 7 and k are the color and flavor indices, respectively) is given by [24J 



A(r,0) 




(12) 



where Ai(r, 6) — \Ai(r)\e ze corresponds to the vortex configuration with winding number one, and A (r) does not 
have a winding number (though it is not spatially constant). If two gaps are the same and constant Aj = A = A^fl, 
it is the CFL phase which is symmetric under the rotation in SU (3) c +l+r- On the other hand, the gap structure 
(fT2")l is invariant only under the transformation in SU (2) c+ l+r <£> E^(i) c +l+r which is a subgroup of S£/(3) c +l+r- 
Still, the symmetry of the bulk CFL phase is restored[37| in the limit r —> 00 because the gap profile function |Ai(r)| 
approaches the value of the CFL phase |Ai(r)| — > |Acfl|- 

For the problem of our interest, it is sufficient to consider the following mean-field Hamiltonian density which 
exhibits the CFL phase: 



i l0 l ■ V - M )V? + [A^(C T C75<) t + h. 



(13) 



An explicit representation of the B-dG equation Hty ■ 



where the gap is now expressed as A^ p oc . . , ■ 



is given by (for a similar representation in the homoge- 



neous CFL phase, see Refs. [2CJ, |3S 
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(14) 



where we have introduced the notation e.g., u r for the quark with color "red" and flavor "up" in the Nambu-Gor'kov 
representation. The matrices T-Lo and A; (i = and 1) are given as follows: 



Ho 
A, 







— 1707- V + fj, ) ' 

Ai7o7 5 
-AJ7075 



(15) 
(16) 
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2. Multiplets in SU(2) C 



+L+R 



Eigenstates of Eq. (fT4|) are classified into multiplets of the unbroken S'?7(2) c+ l+r symmetry: triplet, doublet and 
singlet states. Note first that nine quark states (3 flavors x 3 colors) can be decomposed by the generators of the 
S , C/(3) c+ l+r symmetry of the CFL phase together with a unit matrix [32l |: 




A\ A a 



(17) 



where Xa (A — l,--- ,8) are the Gell-Mann matrices normalized as ^(XaXb) = 2Sab and Ag = y/2/3 • 1. Thus, 
A = 1, • • • ,8 components are the SU(3) octet, and A = 9 the singlet. However, in the presence of a non-Abelian 
vortex, it docs not make sense to mention the multiplets under the S'?7(3) c +l+r transformation. We rather treat 
multiplets under the unbroken S'[/(2) c+ l+r rotation. The triplet and the singlet are respectively given by 

*t = * (1) Ai + * (2) A 2 + * (3) A 3 , (18) 
* s = * (8) A 8 + * (9) A 9 . (19) 

By using the color-flavor representation, each component that appears above is expressed as 'I'W = (d r + u g )/y/2, 
^( 2 ) = (d r -u g )/(V2i), *< 3 ) = (u r -d g )/V2 for the triplet and *W = (u r + d g - 2s b )/V6, *< 9 ' = (u r + d g + s b )/y/3 
for the singlet. In addition to the triplet and singlet states, there are two 'doublet' states. However, we do not discuss 
them here because they do not contain zero modes. On the other hand, as we will see later, the triplet and singlet 
states include zero modes and are thus important. 

One can explicitly check that ^>t and ^f s defined above indeed transform as a triplet and a singlet, respectively. 
The 5?7(2) c+ l+r rotation acts on the quark field ^ as 

* -> = U F * Uj (20) 

where Up = e l9 - x / 2 and U c = e l ^' A / 2 are the 5{/(2) l /r and SU(2) C rotations, respectively. Since we used the vector 
A = (Ai, • • • , As), the parameters 9 and </> are defined only for the first three components 9 — (9x, #2, #3, 0, 0, 0, 0, 0) 
and 4> = ((pi, </>2, </>3, 0, 0, 0, 0, 0). For the SU(2) color-flavor locking, the simplest choice is given by (pi = —9i, (p2 = #2, 
and 03 = —63. Note also that 9i (i = 1, 2, 3) are matrices in the Nambu-Gor'kov representation 



-1 ) ' 62 - ° 2 ( 1 ) ' 03 ~ 03 ( -1 



(21) 



with 9i (i = 1,2,3) being real numbers. Then, one finds that the infinitesimal changes of the triplet ^ t components 
are closed within the three components (i = 1, 2, 3): 

S¥ r) = 9 3 ^ {2) - 6> 2 * (3) , 

s¥ 2) = 0!* (3) - e»3* (1) , (22) 
s¥ 3) = 9 2 ^ {l) - di¥ 2 \ 

and the singlet is invariant 

5^=S^ i9) =0. (23) 



3. Zero-mode solutions 



Apart from the complication coming from the color and flavor degrees of freedom, the structure of the B-dG equation 
is essentially the same as in the case of the single flavor fermion. Namely, requiring fmiteness of the solutions with 
the help of the "Majorana-like" condition, one obtains the zero energy solutions both in the triplet and singlet states. 
The B-dG equation naturally yields both the left and right modes and the structure of the solutions are common in 
both cases as in the case of the single flavor. However, we discuss below only the right-hand modes of the zero modes 
firstly because the chiral symmetry is indirectly broken in the CFL phase (even though it is the SU(2) CFL), and 
secondly because if the system couples to some external fields (as the edge states in the quantum hall effect) only one 
mode (left or right, depending on the properties of the external field) will remain as a zero mode. 



() 




FIG. 1: Distributions of the triplet and singlet zero modes around a non-Abelian vortex. 



The zero modes of the triplet right-handed quarks are analytically given by 

^'m) - c ' ( t:t\ ) • * <2, <--'"> - °' ( -1m) ) ■ * ( "<^> - c > ( ) ■ < 24 > 

where C, are normalization constants and the particle (<£>) and hole (rj) components are 

rfr, 9) = C fo \^\dr' ( JM \ t ^ e) = e - X |A l( ,0|d,' ( ) , (25) 



for generic shape of the vortex profile |Ai(r)|. Notice that these solutions do not contain the unwinding gap |Ao(r)|. 
This is easily understood from the structure of the B-dG equation fH|) . For example, if one looks at the d r and 
u g sector, the Hamiltonian contains only Ai. Another peculiarity about this solution is a minus sign of the hole 
component in \f( 2 ). However, the minus sign must be there so that the zero-mode solutions in the triplet states 
transform as the adjoint representation under the S'J7(2) c+ l+r rotation as shown in Eq. (f22"|) . 

On the other hand, the singlet zero modes depend on both |Ai(r)| and |Ao(r)|, and due to this complication, we 
did not find explicit solutions for generic shape of the gaps. Still, approximate solutions are available where two gaps 
| Ai (r) | and |Ao(r)| are both replaced by the constant gap |Acfl| in the bulk CFL phase: 

where D is a normalization constant and the particle and hole components are given by 

. _ iA CFL i r / J (^r) \ , . _ i^cF L i r / -J 1 (^r)e~' 9 \ 

2 {iJ^e")' m{r > 0) " 6 { J' (27) 

, m - i A cF L i r / J Q (ur)e~ l8 \ , n , _ i^c FL i r / -Ji(fj,r) \ ,„ oN 
^ {r ^ ~ 6 2 { iJ^r) )>V2(r,e)^e * r [u^JiO ) ■ 

These approximate solutions are in fact asymptotic forms of the solutions valid at large distances much away from 
the vortex line. 

In the triplet zero- mode solutions (|25[) . if one approximates the gap profile |Ai| by the constant gap |Acfl| in the 
bulk CFL phase, one finds that the solutions decay as e - i AcFL ' r . On the other hand, the singlet zero modes decay as 
e~l AcPL l r / 2 . Therefore, distribution of the singlet zero modes is wider than that of the triplet, as schematically shown 
in Fig. [Q 



III. LOW ENERGY EFFECTIVE THEORY OF GAPLESS FERMIONS 



In the previous section, we have discussed strictly zero-energy eigenstates of the B-dG equation, and found that 
such zero modes are well-localized on the plane perpendicular to the direction of the vortex line. While the transverse 
motion of such zero modes is frozen, they can move along the vortex line to give a linear dispersion with respect to 
the longitudinal momentum k z . In this section, we develop an effective description of such low-energy excitations. 
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A. Single flavor 



We again start with the single flavor case. We first notice that the Hamiltonian ([2]) can be decomposed into 
"transverse (_L)" and "longitudinal (z)" parts: 

H = f-id x -V±-fi |A|e ie 7o 7 5 \ / -ia z -§- z 
V -|A|e- 4e 7o7 5 ~id ± -V ± +fi) \ 

= n^ + n z . (29) 

Since the longitudinal Hamiltonian H z is linear with respect to the z derivative, it immediately implies that the 
low-energy excitations have linear dispersion relations with respect to k z . Recalling that the zero- mode solutions 
^o{ r , 9) (with k z — 0) are in fact eigenstates of the transverse Hamiltonian (7i± +'H z ) 1 $ l o{r 7 0) = H± i &o(r, 0) — 0, one 
can see that the z dependence enters in a factorized way: 

V(t,z,r,6) = a(t,z)^ o (r,0), (30) 

where we have recovered time dependence and consider the Schrodinger equation id^/dt = T~L^ . Note that the time 
dependence appears only in a(t, z) because the zero-energy states are 'static'. Plugging Eq. (|30[) into the Schrodinger 
equation, and integrating it over transverse coordinates after multiplying it by the zero-mode solution from the 
left, we obtain the equation of motion for a(t, z): 

i—a(t,z) = I *l{r,0)HMt,z)Mr,0)rdrdO, (31) 

where we have used the normalization J ^Q^Qrdrd9 = 1. We can rewrite this equation as 

-+«( M ,|A|)— \a(t,z) = 0, (32) 

with the "velocity" v(n, |A|) defined by 

V ( M ,|A|) = J *S(r,d) ( a Q z ^ ) Mr,9)rdrd6. (33) 

The solution to Eq. (f32|) is given by a(t, z) oc ^ t - lk ^ z with a linear (gapless) dispersion with respect to k z : 

£ = v(fi,\A\)k z . (34) 



Therefore, low-energy excitations inside the vortex are gapless (massless) fermions described by Eqs. (|32|) and (|33l) . 
One can indeed express these fermions in terms of spinors in 1+1 dimensions, and write down an equation similar to 
the Dirac equation. 

One can compute the velocity by using the explicit form of the zero- mode solutions, i.e., the right mode (Eqs. ([5]), 
((6])) and the left mode (Eqs. |7]), (8J). If we approximate the gap profile function |A(r)| by a constant |A|, we obtain 
for the right mode 

2 E( ^ '., ) 

v+( f i,\A\) = — £ 5 ^ — -1, (35) 



and for the left mode 



V-((M,\A\) = -v + (»,\A\) . (36) 



Here, K{x) and E(x) are the complete elliptic integrals of the first and second kinds. Therefore, the right (left) mode 
moves towards the plus (minus) direction of the z axis with the velocity less than the speed of light. In fact, as shown 
in Fig. the velocity v+(/j,, |A|) (|5S")) is a decreasing function of n/\A\, and v + = 1 (the speed of light) in vacuum 
(/i = 0) [the result in vacuum agrees with that obtained by Witten [17| in the context of cosmic strings] . This is also 
seen from the asymptotic forms of the velocity: for small A*/|A| 

« + (/i,|A|)^l-^ + 0(( M /|A|) 4 ), (37) 



s 




FIG. 2: Velocity of gapless fermions along the vortex axis. 



while for large ///A, 

v+fa, |A|) ~ (-1 + 2 In 2 + In J±-^ + 0((|A|/ M ) 4 ) . (38) 

Since the gap in the color superconductivity is parametrically given by |A(//)|//i oc e~ b /s{^) with b being a numerical 
constant and g(n) the running coupling, n/\A\ increases with increasing fi. Therefore, we conclude that the effective 
velocity v+(/j) decreases with increasing /i. If we take typical values /i — 1000 MeV and |A| = 100 MeV, then the 
effective velocity is estimated as v ~ 0.027 which is considerably smaller than the speed of light. 



B. Color-Flavor-Locked phase 

For the zero modes trapped in a non-Abelian vortex, one can perform the same procedure to obtain the effective 
theory. The low-energy effective theories of the left and right fermions in the triplet and singlet states are 

*(|+^)^(M)=0, (39) 

with the dispersion relations 

£ = vfk z , (40) 

where i = t,s are for triplet and singlet, and +, — are for the right and left modes. If we approximate the gap by a 
constant value in the bulk CFL phase, then the velocity of the triplet is exactly the same as the result of the single 
flavor (Eq. (|35|) ) with the replacement of |A| by |Acfl|, ke., 

U t ± = ±U + ( M ,|A C FL|). (41) 

On the other hand, since the singlet zero modes have a wider transverse dependence, the velocity is evaluated by 
using the same function v + in Eq. (|55|) . but with the replacement of |A| by |Acpl|/2, i.e., 

v?=±v + (n,\A CFL \/2). (42) 

In the vacuum /i = 0, both the velocities of the triplet and the singlet are the same as the speed of light. However, 
with increasing densities, the two velocities start to deviate. For example, if one uses the asymptotic form of the 
velocity at large /x/|A| (Eq. (|3"5)0 . then one finds that the velocity of the singlet is about quarter of that of the triplet. 
Therefore, one can draw a schematic picture of the dispersion relations for all the SU(2) multiplets as shown in Fig. [3] 
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FIG. 3: Dispersion relations between energy E and momentum in z direction k z for the ground states of triplet, singlet and 
doublet. 



We studied the fermion structure of a non-Abelian vortex in color superconductivity. By analyzing the Bogoliubov- 
de Gennes equation with the vortex gap profile, we found fermion zero modes at the Fermi level in the center of the 
vortex. We first demonstrated the single flavor case with an Abelian vortex as a warm-up before the more complicated 
case, and then turned to the CFL case with a non-Abelian vortex. We found that there are two gapless modes in the 
triplet and singlet states of S'?7(2) c+ l+r &> f7(l) c +L+R which is the unbroken symmetry in the non-Abelian vortex. 
We also constructed the low-energy effective theory of the gapless fermions which describes the dynamics of fermions 
living on the vortex line (and thus in 1+1 dimensions). The gapless modes propagate along the vortex axis at the 
velocities smaller than the speed of light. 

In this paper, we did not consider the presence of any gauge field in the vortex. However, the ordinary vortex of 
type-II superconductivity can trap a (quantized) magnetic flux in it. In the case of the non-Abelian vortex, there 
exist gauge fields which are associated with the unbroken S'[/(2) c+ l+r (S> t/(l) c +L+R symmetry. However, they are 
massive even in the vortex because only the winding component of gaps vanishes at the center of non-Abelian vortices, 
A(r = 0) = diag(0, *,*) with non-zero constants where all generators of gauge transformations remain broken 
there. Thus they do not affect the low-energy excitations at the energy scale below the masses of gauge bosons. 
Therefore, we expect our effective theory of the gapless fermions is valid at least for very low-energy excitations. 

In the discussion above, we assumed the limit of massless fermions for all flavors. In reality, the strange quark 
has a finite current mass, and cannot be regarded as a massless fermion. The effect of the strange quark mass on 
non-Abelian vortices has been studied recently in the Ginzburg-Landau model (33[. In the B-dG equation for vortices, 
we expect the energy-momentum dispersions of singlet and triplet zero modes are affected differently by this effect 
since the strange quark differently enters them (see Eqs. (fT5)l . (jTHJ) ) - 

In the Ginzburg-Landau approach, a non-Abelian vortex has the CP 2 bosonic zero modes depending on the flux as 
explained in Introduction [27j. A relation between those bosonic zero modes and fermionic zero modes in the B-dG 
equation studied in this paper is unclear at this stage. We expect that it will be clarified if we solve the B-dG equation 
self-consistently, which remains as an important future problem. It has also been studied 27] in the Ginzburg-Landau 
model that there exists a repulsive force between two parallel non-Abelian vortices, and that it does not depend on 
the CP 2 zero modes if they are placed at a distance much larger than their core size. This situation is justified for 
instance in a vortex lattice with the lattice spacing larger than their core size. When a lattice spacing is comparable 
with the core size, we have to use the B-dG equation for multiple vortices. We expect that such an analysis gives a 
force depending on the CP 2 zero modes, namely, on color fluxes which are carried by the vortices. 

Note added: Just before we were about to finish the manuscript, we noticed through correspondences with Y. 
Nishida that he was doing very similar calculations as ours. It turned out from his recent work [34j that his main 
emphasis was on the topological aspects of the fermion zero modes, which we did not discuss in detail. Although his 
analysis is largely the same as ours in the single flavor case, he discussed only the Abelian vortices in the CFL phase, 
which is different from our analysis on the non-Abelian vortices. 



IV. SUMMARY AND DISCUSSION 
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Since the energy of a single vortex A = |A| e'" 9 with winding number |n| > 2 is proportional to n 2 , it decays to n vortices 
each of which has a unit vorticity. Therefore, we consider only the case with n — 1. Note also that the anti- vortex with 
n = — 1 gives qualitatively the same results. 

The Majorana condition with re = 1 leads to the right mode, while re = — 1 the left mode. 

Precisely, the gap approaches A(r — > oo) = AcFLdiag.(e lS , 1, 1), and the extra phase in the first component can be factored 
out by a regular gauge transformation as A — > UA = Acfl e l9 ^ 3 diag.(l, 1, 1). Here U (r, 0) is a regular gauge transformation 
given by U(r, 9) = exp[z(?(r)(#/3)diag.(— 2, 1, 1)] £ 517(3) where g(r) is an arbitrary function with the boundary conditions 
g(r = 0) = and g(r ->■ oo) = 1 



